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I.  INTRODUCTION 


Research  conducted  over  the  past  decade  has  demonstrated  that 
electromagnetic  propulsion  is  a  viable  technology  for  accelerating 
macroscopic  projectiles  to  very  high  velocities.  Particularly 
noteworthy  are  the  achievements  in  rail  gun  technology,  whose  modern 
history  began  with  the  work  of  Marshall  and  his  co-workers  [1-3]  at  the 
Australian  National  University.  The  state-of-the-art  in  rail  gun 
technology  is  exemplified  by  the  recently  completed  program,  sponsored 
jointly  by  DARPA  and  the  U.S.  Army  Armament  Research  and  Development 
Command,  in  which  a  five-meter  rail  gun  built  by  Westinghouse  was  used 
to  accelerate  a  317— g  projectile  to  4.2  km/ s  [4]. 

Although  this  particular  experiment  utilized  a  solid  armature, 
it  is  generally  recognized  that  the  push  to  ever  higher  velocities  will 
require  the  use  of  plasma  armatures,  both  to  circumvent  the  problem  of 
solid  armature  breakup  and  melting  at  high  driving  currents  and  to  help 
insure  good  contact  between  the  armature  and  rails  at  high  projectile 
velocities.  Consequently,  considerable  effort  has  been  expended  in  the 
theoretical  and  experimental  study  of  arc-driven  rail  guns. 

Included  among  the  theoretical  studies  of  arc  armatures  are  the 
analyses  of  Powell  and  Batteh  [5-7],  which  are  based  primarily  on  two 
assumptions  -  (1)  that  the  arc  is  steady  in  a  reference  frame  that 
accelerates  with  the  arc/projectile  system,  and  (2)  that  the  rail  height 
is  mnch  greater  than  the  rail  separation.  The  first  assumption  is 
probably  reasonable  once  the  effects  of  the  initiation  sequence  used  to 
generate  the  arc  have  dissipated.  The  second,  however,  is  decidedly  not 
valid  since  rail  gun  geometries  typically  have  comparable  values  for  the 
rail  spacing  and  rail  height. 


1.  J.P.  Barber,  'The  Acceleration  of  Hacroparticles  and  a  Hypervelocity 
Electromagnetic  Accelerator,'  Ph.D.  Thesis  (Australian  National 
University,  1972)  (Unpublished). 

2.  R.A.  Marshall,  'The  Australian  National  University  Rail  Gun  Project,' 
Atomic  Energy  16.,  January  1975. 

3.  S.C.  Rashleigh  and  R.A.  Marshall,  'Electromagnetic  Acceleration  of 
Macroparticles  to  High  Velocities,'  J.  Appl.  Phys.  ££,  2540  (1978). 

4.  'Laboratory  Demonstration  Electromagnetic  Launcher  -  EMACK, '  Commis¬ 
sioning  Test  Results,  Westinghouse  Research  and  Development  Center, 
Pittsburgh,  PA. 

5.  J.D.  Powell  and  J.H,  Batteh,  'Plasma  Dynamics  of  an  Arc-Driven, 

Electromagnetic  Projectile  Accelerator,'  J.  Appl.  Phys,  2717 

(1981). 

6.  J.H,  Batteh,  'Analysis  of  a  Rail  Gun  Plasma  Accelerator,'  Science 
Applications,  Inc.,  Ballistic  Research  Laboratory  Contract  Report  No. 
ARBRL-CR-00478,  1982  (AD  A114043). 

7.  J.D.  Powell  and  J.H.  Batteh,  'Two-Dimensional  Plasma  Model  for  the 
Arc-Driven  Rail  Gun,'  J.  Appl.  Phys.  5£,  2242  (1983). 


The  two  major  objectives  of  the  effort  conducted  under  this 
contract  were  (1)  to  develop  a  simple  method  for  modifying  the  analyses 
in  Refs.  5-7  to  account  for  finite-height  rails,  and  (2)  to  determine 
the  extent  to  which  including  this  effect  alters  the  predictions  of  arc 
properties.  A  third  objective  was  to  review  the  studies  of  heat 
transfer  and  erosion  in  conventional  gun  tubes  and  assess  their 
applicability  to  the  study  of  damage  to  the  rails  in  rail  guns.  Another 
task  undertaken  during  this  contract  period  was  a  joint  study  with  John 
Powell  of  the  Ballistic  Research  Laboratory  of  the  limitations  imposed 
on  projectile  acceleration  by  the  atmosphere  in  a  non-evacua ted  rail 
gun.  Since  the  results  of  this  latter  task  will  be  published  elsewhere 
181,  only  the  results  of  the  first  three  tasks  will  be  presented  in  this 
report. 


The  report  is  organized  as  follows:  In  Section  II  we 

investigate  the  effect  of  finite-height  rails  on  the  net  accelerating 
force  in  the  rail  gun,  based  on  the  integration  of  the  Lorentz  force 
over  the  cross  section  of  the  armature.  From  this  analysis  we  derive  an 
expression  for  the  effective  inductance  per  unit  length  of  the  gun, 
which  is  a  function  of  the  geometry  of  the  rails  and  armature.  Our 
equation  for  the  inductance  is  actually  more  general  than  Niw&'s  often 
used  formula,  found  in  the  text  by  Grover  [9],  since  we  allow  the  height 
of  the  armature  on  the  rails  to  be  different  from  the  rail  height.  In 
Section  III  we  develop  a  simple  method  for  incorporating  the  effects  of 
finite-height  rails  into  the  analysis  of  rail  gun  arcs.  In  addition,  we 
show  that  for  typical  rail  gun  geometries,  the  momentum  equation  which 
describes  the  pressure  variation  in  the  arc  can  be  well  approximated  by 
merely  scaling  the  corresponding  equation  obtained  for  the  case  of 
infinitely  high  rails.  We  exploit  this  result  in  Section  IV  where  we  use 
a  modification  of  the  analysis  of  Ref.  6  to  determine  the  effect  of 
finite-height  rails  on  arc  properties.  Section  V  contains  a  brief 
discussion  of  the  mechanisms  responsible  for  gun  tube  erosion,  and 
contrasts  these  with  the  mechanisms  most  likely  to  result  in  damage  to 
the  rails  in  arc-driven  rail  guns.  This  section  also  contains  a 
discussion  of  some  recent  models  developed  to  analyze  mechanical  and 
thermal  damage  in  rail  guns.  Section  VI  contains  a  summary  and 

discussion  of  the  work  presented  in  this  report,  as  well  as 
recommendations  for  future  studies. 


8.  J.D.  Powell  and  J.H.  Batteh,  'Atmospheric  Effects  on  Projectile 
Acceleration  in  the  Rail  Gun,'  J.  Appl.  Phys  (in  press). 

9.  F.W.  Grover,  Inductance  Calculations  (Van  Nostrand,  New  Tork,  1946) 
Ch.  10. 
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II.  EFFECTIVE  INDUCTANCE  CALCULATION 


In  this  section,  we  present  a  derivation  for  the  effective 
inductance  per  unit  length  (hereafter  referred  to  as  the  inductance)  of 
a  rail  gun  with  finite-height  rails.  This  inductance,  denoted  by  L,  is 
defined  such  that  the  net  accelerating  Lorentz  force  on  the  armature  and 
projectile  is  given,  for  a  fixed  geometry  and  constant  current,  by 

F  -  LI  *  12  (2-1) 

o 

where  I  is  the  magnitude  of  the  current  in  the  rails, 
o 

The  geometry  used  in  the  calculation  is  shown  in  Fig.  1.  In  the 
figure,  denotes  the  distance  from  the  power  supply  to  the  armature,  w 

is  the  rail  separation,  and  h  is  the  height  of  the  armature  on  the 

a 

rails,  which  in  general  will  be  different  from  the  rail  height,  h  . 

r 

For  typical  rail  gun  configurations,  the  length  of  the  gun  far 

exceeds  the  rail  height  and  separation.  Thus,  for  most  of  the 

acceleration  process  x  >)  h  and  w  and  the  rail s  can  be  assumed  to 

o  r 

extend  to  -  *  in  the  x-direction.  We  incorporate  this  assumption  into 
the  analysis  and  assume,  in  addition,  that  the  current  is  carried  on  the 
inner  surface  of  the  rails.  The  x-component  of  the  current  in  the  rails 
is  taken  to  be  uniformly  distributed  on  the  rails. 

The  armature  is  assumed  to  carry  a  current  only  in  the 
y-direction,  the  current  being  distributed  uniformly  along  the  inner 
surface  of  the  armature.  Thus,  we  can  write  the  armature  current 
density  as 


J  (x,y,z)  =  J  (y,z)  6(x-x  )y  (2-2) 

a  ay  o 

where  &  is  the  delta  function  and 


I  /h  #  lyli  w/2  ,  Izl  <.  h  / 2 
o  a  a 

0  ,  otherwise 

Calculating  the  inductance  based  on  a  current  which  flows  only 

on  the  surface  of  the  armature  is  not  as  restrictive  as  one  might  guess. 

In  Ref.  10  it  was  shown  that  for  h  =  h  =  «*>  the  net  accelerating 

a  r  ° 


10.  J.D.  Powell  and  J.H.  Batteh,  'Plasma  Dynamics  of  the  Arc-Driven  Rail 
Gun,'  Ballistic  Research  Laboratory  Report  No.  ARBRL-TR-02267 , 
1980  (AD  A092345) . 
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pressure  on  the  armature*  and  thus  the  effective  inductance*  are 
independent  of  the  distribution  of  the  armature  current  in  the  axial 
direction  (x-direction) •  A  similar  result*  based  on  the  analysis  of 
Section  HI  of  this  report*  is  derived  in  the  Appendix  for  the  case  where 
hft  and  h^  are  equal  but  finite.  Ve  have  not  attempted  to  prove  this 

result  for  the  more  general  case  where  h  ^  h^*  but  these  analyses 
suggest  its  validity. 

The  current  density  in  the  rails  is 


Jr(x,y,z)  =  [&(y-w/2)  -  6(y+w/2)J  [j^x,  x)x  +  Jrz(x*zVsTj  .  <2—3 ) 

Because  of  our  assumption  that  the  x-component  of  the  current  is 
uniformly  distributed  on  the  rails* 


,  x  <  x  *  1 z|  1  h  12 
o  r 

,  otherwise  .  (2-4) 

J  is  obtained  from  the  current  continuity  equation* 
xz 

"v  •  7-  0  .  (2-5) 

In  writing  Eq.  (2-5),  we  have  neglected  the  contribution  of  the 
displacement  current*  which  is  valid  for  nonrelativistic  projectile 
velocities.  Substituting  the  equations  for  and  into  Eq.  (2-5) 

yields  an  ordinary  differential  equation  for  J  •  Solving  this  equation 

subject  to  the  conditions  that  J  be  continuous  at  z  9  h  /2  and  vanish 
at  z  -  0  yields  rz  a 


<*,*>  = 


I  /h 
o  r 


J  (x*z)  -  -  I/h 
rz  o  a 


6(x-x  ) 
o 


0 

-  h  /h  (z  +  h  /2)* 
a  r  r 

(1  -  \/hr)z 


1  z |  >  h  72 
r 

-h  /2  i  z  i  -  h  /2 
1  a  (2-6) 
Izl  <  h 

a 


ha/hr  (hr/2  -  *) 


.  h.  n  i  z  <  k  /2 

a  r 


Actually*  J  is  presented  here  only  for  completeness  since  it  will  not 

contribute  to  the  accelerating  force  because  of  our  assumption  that  the 
current  in  the  armature  has  only  a  y— component.  Thus*  only  the 
component  of  the  magnetic  induction  in  the  z-direction,  B  *  can 
contribute  to  a  force  in  the  armature  in  the  x-direction. 
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To  calculate  the  effective  inductance,  we  first  calculate  the 
net  force  in  the  z-direction  on  the  armature  from  the  Lorentz  law,  which 
for  our  model  reduces  to 


F 

x 


T  w/2 

I 


h  /2 
a 


B  (x  ,y,z)  dz  dy 
z  o 


-w/2  -h  n 

a 


(2-7) 


and  then  use  Eq.  (2-1)  to  obtain  L.  In  this  report,  we  shall  use  the 
notation  <A(x)>  to  denote  the  average  of  A  over  the  plane  of  the 
armature  located  at  x,  i.e.. 


<A(x)>  = 


wh 


w/2 


I 

-w/2 


ha/2 

I  A(x,y,z)  dz  dy 

-h  12 
a 


(2-8) 


Using  this  notation,  we  can  write  F^  as 


F 

x 


<B  (x  )> 
z  o 


and  the  effective  inductance  as 


(2-9) 


L  =  2w  <B  (x  )>/  I  . 

z  o  o 


(2-10) 


Since  the  armature  current  cannot  exert  a  net  force  upon  itself, 

only  the  rail  currents  need  to  be  considered  in  deriving  <B  (x  )>. 

z  o 

Calculating  B^  due  to  the  rail  currents  for  a  point  on  the  surface  of 
the  armature  from  the  Biot-Savart  law  yields 


B  (x  ,y,z)  = 
z  o  4nh 


h  /2 


-h  12 
r 


f  0  r  *  r  y-w/2 

J  J  L  v 


y+w/2 


dz'dx' 


(2-11) 


—  7 

where  fi  =  4n  x  10  H/m  is  the  permeability  of  free  space,  and  r^  and  r^ 
are  given  by 


14 


x')*  +  (y  -  w/2) 2  + 
x')  +  (y  +  w/2)  + 


(2-12) 


Carrying  ont  the  integration  yields 


*lIo 

Wy>1) "  4^ 

X 

r  -» 

!  tan 

(£■)- 

— i 

+  tan 

/2z+h\ 

--(£■)] 

(2-13) 

\2y+w  /, 

Averaging  Eq.  (2-13) 

over  the 

armature 

cross  section 

and 

substituting 

into  Eq.  (2-10)  yields  the 

following 

expression  for  the  effective 

inductance 

2 

L  = 

2pw  r 

nh  h  1 
r  a  u 

■  fv*- 

A  \  2w 

)  '  1  (  2. 

(2-14) 

where 


I(s)  =  s  tan-1  (s)  +  s*/4  In  (1  +  l/s»)  -  1/4  In  (1  +  s»). (2-15) 

As  mentioned  previously,  Niwa  has  calculated  an  expression  for 
the  effective  inductance  of  a  rectangular  current  sheet,  which  is 
presented  in  the  text  by  Grover  [93.  Our  result,  Eq.  (2-14),  agrees 
numerically  with  Niwa's  formula  when  h^  is  equal  to  hft  and  when  Niwa’s 

formula  is  applied  in  the  limit  of  infinitely  long  rails. 

It  is  apparent  from  Eqs.  (2-14)  and  (2-15)  that  the  inductance 
is  a  function  only  of  the  rail  height  ratio,  h  /w,  and  the  arc  height 

ratio,  kft/hr.  This  functional  dependence  is  depicted  in  Fig.  2  where  we 

have  plotted  the  inductance  as  a  function  of  the  rail  height  ratio  for 
h  /h  <=  0.5  and  1.0.  Also  marked  in  the  figure  are  the  theoretical 

ft  X 

predictions  for  the  geometry  of  the  Rashleigh-Mar shall  (RM)  experiment 

[3]  (h  / w  =  1.5,  h  /h  =  0.67)  and  the  Ballistic  Research  Laboratory 
r  a  r 

(BRL)  experiments  [11)  (h  /w  =  2.0,  h  /h  =  0.5). 

X  ft  X 
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Effective  Inductance,  L  -  u  H/m 


Figure  2  indicates  that  considerable  error  can  be  introduced  if 

Niwa's  formula,  which  corresponds  to  the  solid  curve  in  the  figure,  is 

not  used  judiciously.  For  example,  if  we  apply  Niwa's  formula  to,,  the 

BEL  tests  by  assuming  that  h  and  h  are  equal  to  the  arc  height,  then 

h  r 

we  obtain  a  value  of  L  equal  to  0.628  pH/m,  which  is  considerably 

greater  than  the  value  of  0.444  pH/m  predicted  for  h  /h  =  0.5. 

ft  r 

Somewhat  better  agreement  is  obtained  if  both  heights  are  assumed  equal 
to  the  rail  height.  But  in  this  case,  one  must  be  careful  to  apply  the 
force  to  the  actual  armature  cross-sectional  area  in  order  to  obtain 
accurate  predictions  for  the  pressure  on  the  armature. 

In  previous  studies,  Powell  and  Batteh  have  attempted  to 

simulate  the  RM  experiment  based  on  a  rail  gun  model  which  assumes,  for 
computational  simplicity,  that  the  rails  are  infinitely  high  [5].  As 
discussed  in  Ref.  6,  this  assumption,  for  a  current  per  unit  height 
equal  to  that  of  the  EM  experiment,  overestimates  the  arc  pressure  and 
acceleration  by  approximately  8  factor  of  two.  In  fact,  to  obtain  the 
same  pressure  as  that  predicted  from  the  inf inite-height  model  from  Eq. 
(2-1)  and  the  actual  RM  geometry  would  require  an  effective  inductance 
on  the  order  of  lpH/m.  Equation  2-14,  on  the  other  hand,  predicts  a 
considerably  lower  value,  0.526  pH/m,  for  the  RM  geometry.  Thus,  we 

expect  that  the  pressures  and  accelerations  predicted  in  the  earlier 
studies  are  inaccurate.  This  point  will  be  addressed  in  greater  detail 
in  the  following  section.  It  is  perhaps  worth  mentioning  that  in  Ref. 
6,  an  effective  inductance  of  0.53  pH/m  was  estimated  for  the  RM 
experiment  based  on  the  experimentally  measured  current  profile  and 
projectile  acceleration.  This  value  is  essentially  identical  to  the 
theoretical  calculation  based  on  Eq.  (2-14).  Ve  suspect,  however,  that 
the  remarkable  agreement  may  be  fortuitous  to  some  extent  since  the 
model  leading  to  Eq.  (2-14)  does  not  account  for  frictional  drag  on  the 
projectile  and  for  the  modification  of  the  inductance  due  to  the 

diffusion  of  current  in  the  rails,  both  of  which  are  implicitly  included 

in  the  calculation  based  on  the  experimentally  measured  parameters. 

The  theory  developed  in  this  section  also  provides  a  measure  of 
the  nonuniformity  of  the  accelerating  force,  which  is  often  an  important 
factor  in  projectile/armature  design.  Because  of  our  assumption  that 
the  current  is  uniform  on  the  surface  of  the  armature,  the  ratio  of  the 
pressure  at  any  point  on  the  armature  surface  to  the  average  pressure  is 
proportional  to  the  ratio  of  the  magnetic  induction  at  that  point  to  the 
average  induction,  i.e., 


p(x  ,y,z)/<p(x  )>  =  B  (x  , y, x) / <B  (x  )>. 
o  o  x  o  z  o 


(2-16) 
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It  can  be  readily  shown  by  differentiation  of  Eq.  (2-13)  that  the  peak 
pressure  occurs  at  z  -  0,  y  -  +  w/2,  Consequently,  the  peak,  to  average 
pressure  on  the  projectile  can  be  written  as: 


p 

max 

p 

avg 


.-HE  — 

2irh  L 
r 


[ 


n  +  2  tan 


(2-17) 


Again,  it  is  clear  that  the  pressure  ratio  depends  only  on  the  rail 
height  ratio  and  the  arc  height  ratio. 


Figure  3  shows  plots  of  the  pressure  ratio,  as  determined  from 
Eq.  (2-17),  which  correspond  to  the  inductance  plots  of  Fig.  2.  Figures 
2  and  3  clearly  demonstrate  the  well-known  fact  that  the  effective 
inductance,  and  thus  the  accelerating  force  for  a  given  current, 
increase  with  decreasing  rail  height  ratio,  but  that  this  increase  is 
obtained  at  the  expense  of  a  greater  nonuniformity  in  the  pressure  on 
the  armature.  On  the  other  hand,  the  figures  indicate  that  for  a  given 
rail  height  ratio,  the  effective  inductance  actually  increases  somewhat 
as  h&/hr  is  reduced,  while  a  substantial  improvement  in  pressure 

uniformity  is  achieved.  While  this  suggests  that  there  may  be  an 
advantage  to  designs  employing  a  low  h  /h  ,  it  should  be  noted  that 

tt  i 

decreasing  h  / h  for  a  given  h  will  also  cause  an  increase  in  p 

a  r  °  r  ^avg 

since  a  larger  force  will  be  applied  to  a  smaller  cross-sectional  area. 

In  addition,  there  will  be  increased  ohmic  losses  per  unit  volume  in  the 

armature  and  rail  as  the  current  constricts  to  adjust  to  the  armature 

height. 


18 


Pressure  Ratio, 


Figure  3.  Pressure  Ratio  for  Finite-Height  Rail  Guns 
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III.  MODELING  THE  ARC  DYNAMICS  OF  FINITE-JIEIGHT  ARC-DRIVEN  RAIL  GUNS 


Previous  analyses  of  arc-driven  rail  guns  have  invoked  the 
assumption  that  the  rail  height,  h^,  and  arc  height,  h&,  are  infinite 

[5-7].  This  assumption  leads  to  a  considerable  simplification  in  the 
momentum  equation  describing  the  arc. 

Consider,  for  example,  the  model  shown  in  Fig.  1  where  we  take 

the  armature  to  be  an  arc  with  length  £  =  x,  -  x  .  We  assume  that  the 

a  l  o 

arc  is  steady  in  a  frame  accelerating  with  the  arc/projectile  system  and 
that  the  thermodynamic  properties  of  the  arc  are  functions  only  of  x. 
Then,  if  the  arc  and  rail  heights  are  infinite,  the  magnetic  induction 
field  is  independent  of  y  and  z  and  the  axial  momentum  equation  in  the 
accelerating  frame  becomes  [5] 

dp/d£  =  £a  [Jy($)Bz($)  -  p($)a]  (3-1) 

where  p  is  the  pressure  in  the  arc,  Jy  is  the  current  density  in  the  arc 

in  the  y-direction,  p  is  the  density,  a  is  the  common  arc/projectile 
acceleration  and 


5  = 

(x  -  xo)/t> 

(3-2) 

is  the 

non-dimensional  distance 

measured  from  the  back 

of  the  arc. 

Furthermore,  the  current  density 

and  magnetic  induction  are 

defined  by 

Jy(t) 

«<«>j0 

£  a 

(3-3) 

- 

a 

X 

Jt  y«» d* 

(3-4) 

where  a 

is  the  conductivity  of  the  arc,  a  is  the  average  conductivity  of 

the  arc 

a  - 

i 

f  at 

J  0 

(3-5) 

and 

^  o 

=  I  /h 
o  a 

(3-6) 
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is  the  current  per  unit  height  of  the  rail 


As  formulated  in  Eq.  (3-1),  the  pressure  for  the  infinite-height 
rail  gun  is  assumed  to  be  a  function  only  of  the  axial  coordinate.  This 
assumption  is  not  valid  for  rail  guns  with  finite-height  since  both  Jy 

and  B_  depend  on  z  and  y,  as  well  as  x.  On  the  other  hand,  Eq.  (3-1) 
z 


will  be  a  valid  representation  for  the  average  pressure  on  a 
cross  section  of  the  arc  in  a  finite-height  gun,  provided  that  the  term 
J  B  in  that  equation  is  replaced  by  the  average  of  the  product  over  the 


y  z 

cross  section.* 
would  entail  a 
dimensions. 


Unfortunately,  the  exact  calculation  of  this  average 
detailed  solution  of  Maxwell's  equations  in  three 


To  avoid  this  complication,  we  proceed  in  the  following  manner. 
First,  we  assume  that  the  current  in  the  arc  flows  only  in  the 
y-direction  and  is  uniform  across  the  arc  cross  section.  We  take  its 
value  to  be  the  average  value  of  Jy  over  that  cross  section.  We  then 

estimate  the  average  pressure  on  the  arc  cross  section  according  to  the 
equation 


d<p>/d$  =  £  [  <J>  <B  >  -  p(()a  3  (3-7) 

a  y  z 

where  <B  >  is  calculated  based  on  the  current  distribution  <J  >. 
z  y 

The  average  arc  current  is  given  by 

<Jy>  *  c(£)  (<Ey>  -  v<Bx>)  ,  (3-8) 

where  v  is  the  material  velocity  of  the  arc  and  projectile  in  the 
x-direction,  and  <Ey>  is  the  average  electric  field  in  the  y-direction. 

<Ey>  is  derived  from  Maxwell's  equation 

TxT  =  -  0B?dt  .  (3-9) 

Recalling  that  our  system  is  steady  in  the  accelerating  fram^  and 
transforming  to  that  frame,  we  find  that  the  z-component  of  Eq.  (3-9) 
reduces  to 

dE  /0$  -  8  dE  /dy  -  v  dB  /d$  .  (3-10) 

y  a  x  z 

It  should  be  noted  that,  although  we  have  transformed  to  the 
accelerating  frame,  the  electromagnetic  fields  in  Eq.  (3-10)  are  those 
measured  by  an  observer  in  the  fixed  frame. 


Actually,  an  additional  term  should  be  included  in  Eq.  (3-1)  consisting 

of  the  product  of  J  and  B  ,  but  this  term  should  be  much  less  than  the 

z  y 

JyB^  product  for  typical  rail  gun  geometries. 
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Averaging 
that  E 
yields  X 


Eq.  (3-10)  over  the  arc  surface  at  each  £-plane  and  noting 
vanishes  at  the  surface  of  the  perfectly  conducting  rails. 


d<E  >/d$  =  v  d<B  >/d$ 

y  * 

The  solution  to  Eq,  (3-11)  is 

<E  >  *  v  <B  >  +  E 

y  z  o 


(3-11) 


(3-12) 


where  E^  is  a  constant.  Substituting  Eq.  (3-12)  into  Eq.  (3-8)  and 
using  the  current  conservation  relation 


yields 


and 


J.  <v«- 


E 

o 


jp  ht 
Ja  ha  5 


<J 


a  _  _  j0  hr  o($)  . 

y  e  _ 


£  h  a 
a  a 


(3-13) 


(3-14) 


(3-15) 


Except  for  the  factor  h  /h  ,  Eq.  (3-15)  is  identical  to  the  result 

r  a 

obtained  for  the  arc  current  density  in  an  infinite-height  rail  gun  [5]. 


We  now  proceed  to  the  calculation  of  <B  >  for  the  current  profile 

z 

<J  >.  For  simplicity,  we  limit  our  discussion  to  the  case  where  h  =  h 
y  r  a 

and  denote  this  common  value  by  h.  Figure  4  shows  a  circuit  which 

encompasses  a  portion  of  the  arc.  The  circuit  is  parallel  to  the  rails, 

has  a  height  corresponding  to  the  arc  height  h^,  and  extends  to  infinity 

in  the  x-direction.  Applying  Ampere's  law  to  this  circuit  and  noting 
that  the  electromagnetic  fields  must  vanish  at  infinity  yields 
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Figure  4.  Circuit  Geometry  for  Magnetic  Induction  Calculation 
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h/2 

-  xx 

J  B  (1,7,1)  dz  a  jih  f  <J  >  dz 

-h/2  z  Jx  7 

U~X6) 

-  2  f  B  (x',y,h/2)  dx*  . 

Jx 

In  writing  Eq.  (3-16)  we  have  made  use  of  the  antisymmetry  in  about 

z-0.  Integrating  Eq.  (3-16)  over  y  and  making  the  change  of  variables 
denoted  by  Eq.  (3-2)  yields 


<v*» 


<J  ($')>  d$' 
y 


2£  r 

hH  J 


w/2 


— w/2 


($\y,h/2)  d^'  dy 


(3-17) 


for  the  average  B  field  on  the  arc  at  fc. 
z 

Since  the  current  in  the  rails  is  in  the  x-direction,  it  cannot 
contribute  to  B^.  Furthermore*  we  again  assume  that  the  rails  are 

infinitely  long  in  the  negative  x-direction.  Then  B^  arises  only  from 

the  current  in  the  arc*  and  is  given  by  the  Biot-Savart  law: 


B^(5*y*h/2)  ■ 


x 


0 


h/2  w/2 


J  J 

-h/2  -w/2 


(3-18) 

<J  <$')>  (h/2-z')  dy '  dz'  d{' 

- xJ! - 

»/a 

[«,(?-{')  +  <y-y')  ♦  (h/2-i')  J 


Carrying  out  the  integrations  over  y'  and  z'  in  Eq.  (3-18)  yields 


Bxt{.y.h/2)  -  J  [f<*/2+y.O) 

+  f(w/2-y.O)  -  £(w/2+y,h)  -  £(w/2-y,h)]  d{* 


(3-19) 


where 
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xl 1 

r  *  *  * 

f(y,z)  =  In  y  +  ($-£')  +  y  +  z  J 

-  1/2  In  +  **] 


(3-20) 


Further  integration  of  Eq.  (3-19)  is,  of  course,  not  possible  without 
specifying  the  current  distribution  in  the  arc. 

Substituting  Eq.  (3-19)  into  Eq.  (3-17)  and  performing  some 
rather  tedious  integrations  yields 


a 

(?)>  =  ,2  f  <J  ($’)>  d?' 

*  y 

i 

-  2SF I.  <y«'»  tt  V*]  «' 


(3-21) 


where 


S(®>  -  S(ti) 


.  A  1  0 


s(«>  -  2s(0)  +  s(-n),  n  <  o 


and  the  function  S(q)  is 


(3-22) 


(3-23) 


S<n)  ■  n  In 


(q*.+  h'») 


1  +  ^1  + 

i  +Vi  +  k'»  +  H* 


+  In  (i\  +\  1  +  na) 

-  (l  -  h'»)  In  (n  +  \  l  +  h'*  +  n*) 

-  h'»  In  (n  +  V h'a  +  i*> 

+  2h'  tan“1(n/h») 

-i  [i  +  W i  +  v*  +  n’  1 
+  2h*  tan  [  - J 


(3-24) 


n  (n  -  \i  +  T!a  -  Vh#*  +  n*  +  V1  +  h#*  +  n*>  . 
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In  the  limit  that  q  goes  to  infinity  and  zero*  S (t| )  takes  on  the  values 


S(»)  =  nh#  +  2h*  tan  (1/h1) 

S(0)  -  nh»  -  (1  -  h'»)  In  1  +  V*  -  h'*  In  h' 


(3-25) 


Equations  (3-7),  (3-15),  and  (3-21)  represent  the  equations  for 
the  average  pressure  in  finite-height  rail  guns.  As  such,  they  would 
replace  Eqs.  (3-1),  (3-3)  and  (3-4)  in  the  analyses  of  Refs.  5  and  6. 
It  is  a  straightforward  exercise  to  show  that  the  two  sets  of  equations 
become  identical  in  the  limit  that  h^  s  h^  s  h  approaches  infinity. 

One  result  of  the  form  of  Eqs.  (3-15)  and  (3-21)  is  that  the  net 

force  on  the  arc/projectile  system,  which  is  proportional  to  the 

integral  of  <J  >  <B  >,  is  independent  of  the  current  distribution  in  the 
y  z 

arc.  (This  characteristic,  which  had  previously  been  noted  for 
infinitely  high  rails  [10],  is  derived  in  the  Appendix.)  In  particular, 
the  net  force  is  the  same  as  if  the  current  was  confined  to  the  surface 
of  the  arc,  which  is  precisely  the  assumption  used  to  derive  the 
effective  inductance  in  Section  II.  Thus,  the  force  can  be  obtained  by 
calculating  L  from  Eq.  (2-14)  and  substituting  the  result  into  Eq. 
(2-1). 


It  is  interesting  to  speculate  whether  the  magnetic  pressure 

profile  for  finite-height  rail  guns  can  be  adequately  approximated  by 

merely  scaling  the  pressure  profile  obtained  for  the  case  of  infinitely 

high  rails.  The  rationale  for  attempting  to  do  this  is  the  following. 

In  arc  models  designed  for  computational  efficiency,  such  as  that 

described  in  Ref.  6,  as  well  as  in  models  for  studying  arc  instabilities 

it  is  often  desirable  to  choose  a  functional  form  for  the  electrical 

conductivity  which  simplifies  the  calculation  of  the  magnetic  pressure. 

It  is  apparent  from  the  relatively  simple  form  of  Eqs.  (3-3)  and  (3-4) 

that  this  can  be  readily  achieved  for  the  case  of  rails  with 

infinite-height.  On  the  other  hand,  one  would  be  hard  pressed  to  find  a 

conductivity  which  would  lead  to  a  simplification  in  evaluating  <B  > 

z 

from  Eq.  (3-21).  If  the  pressure  profiles  do  scale,  however,  then  the 
analysis  for  finite-height  rails  would  be  a  straightforward  extension  of 
the  analysis  of  infinite-height  rail  guns. 

From  Eq.  (2-1)  and  the  definition  of  j  ,  we  can  write  the  average 
pressure  over  the  arc  due  to  the  electromagnetic  forces,  which  we  denote 
as  the  magnetic  pressure,  as 
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(3-26) 


In  the  limit  that  h  =  h  *=  *, 

a  r 


L  j 


*  h  * 
o  r 


2  h  w 
a 

reduces  to 


(3-27) 


which  suggests  a  scaling  factor,  f ,  given  by 

-  -  1  V 

f  =  “  “"T"w  (3-28) 

ffl  BP  u  h  w 

.  a 

for  the  same  The  question  then  becomes  how  well  the  magnetic 

pressure  gradient 


d<p>/d£  =  £  <J  ><B  > 

m  a  y  z 


(3-29) 


is  approximated  by  the  scaled  equation 

d  P-s'**  =  'a  f  Jy  Bz  (3-30) 

where  J  and  B  are  solutions  of  Eqs.  (3-3)  and  (3-4). 
y  z 

To  address  this  question,  we  consider  the  case_  where  the 

conductivity  is  constant  throughout  the  arc  and  equal  to  a.  Ve  also 

restrict  our  discussion  to  equal  rail  and  arc  heights  since  that 

restriction  was  imposed  in  deriving  Eq.  (3-21).  For  that  case. 

Equations  (3-3)  and  (3-4)  predict  that  J  is  constant  and  B  varies 

y  z 

linearly  with  g.  If  we  substitute  the  solutions  into  Eq.  (3-30)  and 
define  a  non-dimensional  magnetic  pressure  according  to 


P 

m 


Pm  W 
in 


L  h  j  : 
a  J  o 


Eq.  (3-30)  for  the  scaled  pressure  becomes 

d  p  /jg  -  i  -  e  . 

ms 


(3-31) 


(3-32) 


This  linear  variation  is  shown  as  the  solid  line  in  Fig.  5.  Also  shown 

in  the  figure  are  plots  of  the  normalized  pressure  gradient  for  two 

finite-height  rail  gun  geometries  corresponding  to  h*  equal  to  unity  and 

8  /w  =  2  and  10.  These  curves  were  obtained  by  solving  Eq.  (3-21) 
a 

numerically  and  substituting  the  result,  along  with  <J  >  from  Eq. 
(3-15),  directly  into  Eq,  (3-29).  ^ 
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Normalized  Magnetic  Pressure  Gradient,  d<P  >/d£ 


Comparison  of  the  corves  in  Fig.  5  indicates  that  the  scaled 
result  appears  to  provide  a  good  representation  of  the  magnetic  pressure 
gradient  in  finite-height  rail  guns  for  typical  geometries  (hf 
approximately  equal  to  1.0,  fi  /w  >  5.0),  at  least  for  the  case  of 

ft  ** 

uniform  conductivity.  In  general,  the  agreement  is  better  for  larger 
values  of  fi  / w.  Although  it  is  not  shown  in  the  figure,  agreement  also 

ft 

improves  as  h’  increases,  which  is  to  be  expected  since  the  scaled  curve 
corresponds  to  the  case  where  h'  =  00 . 

Figure  6  shows  the  same  set  of  curves  for  the  case  where  the 
current  profile  in  the  arc  is  assumed  to  vary  linearly  according  to  the 
relation 


Jy  -  -  2j0  ?/Ia  .  (3-33) 

The  current  variation  described  in  Eq.  (3-33)  is  not  chosen  to 
correspond  to  any  realistic  physical  situation  but  merely  to  show  how 
the  current  profile  affects  the  accuracy  of  the  scaled  pressure 
equation.  Again,  we  note  that  the  scaled  pressure  gradient  profile 
gives  a  good  approximation  to  the  numerically  calculated  profile, 
particularly  for  large  values  of  fi  / w. 

ft 

An  interesting  feature  of  the  exact  curves  in  Figs.  5  and  6  is 
that  the  magnetic  pressure  gradient  actually  becomes  negative  near  the 
projectile,  corresponding  to  a  reversal  in  the  direction  of  the  B  field. 
This  arises  because,  near  the  leading  edge  of  the  arc,  the  contribution 
to  the  pressure  gradient  from  the  current  in  the  rails  is  exceeded  by 
the  contribution  from  the  arc's  own  current,  which  tends  to  'pinch'  the 
arc  in  the  x-direction. 
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Normalized  Magnetic  Pressure  Gradient,  d<P  >/d£ 


IV.  EFFECT  OF  FINITE-HEIGHT  RAILS  ON  ARC  PROPERTIES 


In  this  section,  we  investigate  the  extent  to  which  including  the 
effect  of  the  finite  height  of  the  rails  alters  the  predictions  of  the 
properties  of  the  arc  in  arc-driven  rail  guns.  To  do  this  we  exploit 
the  result  of  the  analysis  in  Sec.  3.  namely,  that  the  magnetic  pressure 
gradient  for  finite-height  rails  can  be  approximated  by  scaling  the 
pressure  gradient  for  infinite-height  rails. 

The  model  we  use  is  the  simple  model  for  arc  properties  derived 
in  Ref.  6.  In  that  model,  the  momentum  equation  appropriate  for 
infinite-height  rails  is  used  to  calculate  the  pressure  profile  in  an 
arc  with  constant  electrical  conductivity  and  constant  temperature.  The 
temperature,  in  turn,  is  taken  to  be  the  average  temperature  in  the  arc 
based  on  the  solution  in  three  dimensions  of  the  radiation  heat  transfer 
equation  with  a  source  term  to  account  for  ohmic  dissipation  in  the  arc. 
Several  other  simplifying  assumptions  are  made  and  these  are  discussed 
in  detail  in  Ref.  6. 

Based  on  the  scaling  arguments  discussed  in  Sec.  III.  we  can 
utilize  the  model  of  Ref.  6  to  study  the  properties  of  finite-height 
rail  guns  provided  that  the  following  substitutions  are  made 

D  — h 

a 

j  ►  j  h  /h  (4-1) 

o  J  o  r  a 

p  — ►L(h  /w) 

& 


wherever  the  parameters  D,  j q  and  p  occur  in  the  model.  The  first  two 

substitutions  are  required  to  obtain  the  correct  ohmic  dissipation  in  an 

arc  where  h  /  h  .  (The  theory  of  Ref.  6  was  based  on  h  —  h  =  D.) 
a  r  r  a 

The  substitution  for  p  assures  that  the  correct  scaling  of  the  magnetic 

pressure,  as  defined  by  Eq.  (3-28),  is  achieved. 


The  specific  conditions  we  analyze  are  those  listed  in  Table  1 
which  are  based  on  the  EM  experiment  [3].  The  parameter  D,  equal  to 

'  "  "  current  per  unit  height  for  the 


(h  h  ) 
r  a 


is  used  to  determine  the 


infinite-height  rail  gun  calculation.  The  model  requires  the 
specification  of  either  the  arc  mass  or  arc  length.  We  have  chosen  to 
use  the  arc  length  since  it  is  the  more  readily  measured  of  the  two 


parameters.  The  estimated  value  of  10  cm  for 
previous  studies  of  the  RM  experiment. 


is  consistent  with 
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Table  1.  CONDITIONS  FOR  HIE  RM  EXPERIMENT 


SYMBOL 

QUANTITY 

VALUE 

w 

Rail  Separation 

1.27  i  10"2m 

h 

a 

Arc  Height  on  Rails 

1.27  x  lO"1* 

h 

r 

Rail  Height 

1.91  x  10”*m 

D 

Effective  Rail  Height* 

_2 

1.56  x  10  m 

m 

P 

Projectile  Mass 

3.0  x  10~*kg 

I 

0 

Current 

300  LA 

m 

0 

Mass  of  Ions  and  Neutrals 

1.1  x  10  *kg 

* 

a 

Arc  Length 

0.10  m* 

1/2 

D  =  (h  h  )  »  used  in  the  infinite  rail  height  model 

a  r 

*  Estimated  value 
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Table  2  shows  a  comparison  of  the  predictions  of  the 
infinite-height  model  and  the  finite-height  model  for  various  arc 
properties.  The  final  values  for  the  arc  pressure  and  density  denote  the 
values  at  the  projectile  surface.  As  expected,  the  finite-height  rail 
gun  model  predicts  significantly  lower  arc  pressures  and  accelerations. 
Since  the  two  models  predict  similar  temperatures,  the  density  predicted 
by  the  finite-height  model  is  also  much  lower  than  that  for  the 
infinite-height  model,  resulting  in  a  lower  arc  mass  for  a  given  arc 
length.  On  the  other  hand,  the  muzzle  voltage  predicted  by  the 
finite-height  model  is  approximately  fifty  percent  higher  than  that 
predicted  by  the  infinite-height  model.  This  increase  is  due  to  the 
smaller  arc  cross-sectional  area  (by  the  ratio  of  h  /D)  used  for  the 

finite-height  calculation,  and  to  a  decrease  of  about  fifteen  percent  in 
the  arc  conductivity. 

Since  there  is  a  large  degree  of  uncertainty  in  the  experimental 

value  of  £  ,  Figs.  7-10  show  graphical  comparisons  of  the  predictions 
& 

of  the  two  models  for  the  arc  pressure,  temperature,  muzzle  voltage,  and 
arc  mass.  These  figures  clearly  demonstrate  that  including  the  effect 
of  finite-height  rails  does  not  alter  the  trends  predicted  by  the  model: 
namely,  that  an  increase  in  arc  length  (or,  equivalently,  arc  mass) 
results  in  decreases  in  the  arc  pressure  and  temperature  and  in  the 
muzzle  voltage  Of  the  rail  gun. 
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Table  2,  COMPARISON  OF  MODEL  PREDICTIONS  FOR  RM  PARAMETERS 


SYMBOL 

Pf 

a 

P 

Pf 

ID 

a 

“e 
V 

o 

T 

T 

max 


QUANTITY 

Average  Arc  Pressure  (Nt/ms) 

Final  Arc  Pressure  (Nt/m*) 
Arc-projectile  Acceleration  (m/ss) 
Average  Arc  Density  (kg/m1) 

Final  Arc  Density  (kg/m1) 

Arc  Mass  (kg) 

-i 

Average  Electron  Number  Density  (m  ) 
Mnzxle  Voltage  (V) 

Average  Arc  Temperature  ( °K) 

Maximum  Arc  Temperature  (°K) 


VALUE 
(Ww  =  «) 

1,48x10s 

2,13x10s 

1,40x10 7 

14.3 

20.6 

0.29x10”* 

1.90x10s* 

43. 

33,500 

53,400 


VALUE 

ly*  °  1,0) 

0.94x10s 

1.39x10s 

0.75xl07 

9.67 

14.3 

0.16x10”* 

1.27x10s* 

63. 

31.800 

50,700 


a 


Average  Degree  of  Ionization 


1.46 


1.44 


-  h  /w  =  1.0 

a 

-  ha/w  =  “ 


Figure  9. 


Comparison 


of  r 
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V.  RAIL  DAMAGE  IN  ARC-DRIVEN  RAIL  GUNS 


The  successful  development  of  a  practical*  repetitively  fired 
rail  gun  requires  an  understanding  of  the  mechanisms  responsible  for 
damage  to  the  rails  and  the  development  of  techniques  to  minimize  this 
damage.  In  this  section  we  discuss  some  probable  causes  of  rail  damage 
and  compare  these  with  the  mechanisms  responsible  for  gun-tube  erosion 
in  conventional  guns.  Some  recent  models  developed  to  analyze  erosion 
in  rail  guns  will  also  be  discussed. 

Although  a  general  theory  has  yet  to  be  formulated*  erosion  in 
conventional  gun  tubes  is  believed  to  arise  from  a  combination  of  the 
mechanical  interaction  of  the  projectile  with  the  gun  tube*  and  the 
heating  of  the  tube  by  the  hot  propellant  gases  [11*12].  Mechanical 
erosion*  which  is  relatively  independent  of  the  gas  temperature*  results 
in  a  small  rate  of  wear  and  is  most  important  for  low-velocity  weapons. 
For  large-caliber*  high  velocity  weapons  thermal  effects  are  generally 
the  major  causes  of  barrel  erosion.  The  cyclic  heating  and  cooling  of 
the  steel  gun  tube  in  a  repetitively  fired  gun  leads  to  a  chemical 
transformation  at  the  surface  which  makes  it  particularly  susceptible  to 
cracking.  If  the  surface  of  the  gun  tube  reaches  a  sufficiently  high 
temperature*  then  the  gun  surface  will  recede  as  it  melts  and  is  wiped 
away  by  the  flowing  projectile  gases. 

It  is  likely  that  thermal  and  mechanical  interactions  will  al so 
be  the  major  contributors  to  rail  erosion  in  rail  guns.  Buckingham  has 
presented  a  study  of  mechanical  erosion  in  arc-driven  rail  guns  [13]. 
His  study  focuses  primarily  on  the  drag  and  mass  loss  of  the  projectile 
as  it  interacts  with  the  rails .  The  calculations  are  based  on 

gram-sized  projectiles  which  are  stressed  beyond  their  yield  values  and 
subjected  to  accelerations  on  the  order  of  107  m/s1.  For  these 
conditions,  he  finds  that  a  projectile  with  a  steel  sabot  will 
experience  significant  mass  loss  due  to  viscous  effects  if  one  attempts 
to  accelerate  it  to  much  beyond  4  km/s.  Substantially  better 
performance  is  predicted  for  a  projectile  with  either  a  graphite  or 
Teflon  sabot*  or  for  a  configuration  where  there  is  initially  a  gap 
between  the  projectile  and  rails.  An  effect  which  was  not  considered, 
and  which  may  be  significant,  is  the  leakage  of  plasma  through  the  gap 
which  could  lead  to  a  current-conducting  path  ahead  of  the  projectile. 
To  seal  the  gaps  and  to  prevent  the  projectile  from  'bouncing1  against 
the  rails#  Buckingham  suggests  the  use  of  ablative  bands  around  the 
proj  ectile. 


11.  Interior  Ballistics  of  Guns.  AMCP  706-150  (US  Army  Materiel  Command* 
Washington*  D.C.*  1965)  Ch.  3. 

12.  J.  Corner,  The  Theory  of  the  Interior  Ballistics  of  Guns  (Wiley,  New 
York,  1950)  Ch.  10. 

13.  A.C.  Buckingham*  'Electromagnetic  Propulsion:  Drag  and  Erosion 

Modeling,'  AIAA  J.  19,  1422  (1981). 


39 


Analogous  to  the  propellant  gases  in  a  conventional  gun,  the 
plasma  arc  is  the  major  contributor  to  rail  heating  in  arc-driven  rail 
guns.  There  are,  however,  several  important  differences.  Firsts  the 
average  arc  temperature  is  typically  20,000  to  40,000  °K,  whereas  the 
temperature  of  the  propellant  gases  rarely  exceeds  a  few  thousand 
degrees.  Because  of  the  higher  temperature,  energy  is  transferred  to 
the  rails  primarily  by  radiation.  Furthermore,  most  of  the  transfer  at 
any  point  on  the  rails  occurs  while  the  arc  is  in  contact  with  that 
point.  A  secondary  source  of  rail  heating  is  ohmic  dissipation  in  the 
finite  conductivity  rails. 

Powell  has  recently  completed  a  study  of  rail  damage  arising  from 
these  two  effects  which  is  based  on  solutions  of  the  time-dependent 
diffusion  equation  [14].  It  is  assumed  in  the  analysis  that  once  a 
layer  of  the  rail  reaches  the  melting  temperature,  it  is  removed  and  can 
no  longer  absorb  radiation.  (This  is  equivalent  to  the  'mel t— and-wipe' 
hypothesis  frequently  used  to  study  gun— tube  erosion.)  The  theory 
provides  predictions  for  the  time  to  initiate  melting  at  the  surface 
and,  once  the  melting  temperature  has  been  attained,  for  the  rate  of 
recession  of  the  rail  surface  as  a  function  of  axial  position  on  the 
rails  and  the  arc  properties.  The  following  discussion  will  help 
identify  the  arc  properties  which  are  most  influential  in  determining 
the  energy  transfer  to  the  rails. 

If  the  rail  thickness  is  large  compared  to  the  electrical  skin 
depth,  the  temperature  rise  at  the  surface  of  the  rails  during  contact 
with  the  arc  can  be  approximated  by  [15] 


AT  =  2qg  [t/(n  pr  Cr  kr)]1/*  +  ATR 


(5-1) 


where 


AT, 


R 


P 


*  P, 


I  * 
o 


C  h  * 
r  a 


(5-2) 


represents  the  temperature  rise  due  to  resistive  heating.  In  Eqs.  (5-1) 
and  (5-2),  is  the  value  of  the  radiation  flux  incident  on  the  surface 

of  the  rails,  t  is  the  time  measured  from  the  arrival  of  the  leading 
edge  of  the  arc  at  the  point  in  question,  and  p^,  C£,  k^  and  a  denote 

the  density,  specific  heat,  thermal  conductivity  and  electrical 
conductivity,  respectively,  of  the  rails.  Equation  (5-1)  is  valid  until 
the  surface  temperature  attains  the  melting  temperature.  In  writing 


14.  J.D,  Powell,  'Thermal  Energy  Transport  from  Arc  to  Rails  in  an 
Arc-Driven  Rail  Gun,'  Ballistic  Research  Laboratory  Report  No. 
ARBRL-TR-02530 ,  Oct  83. 

15.  J.D.  Powell,  private  communication. 
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this  equation,  we  have  noted  that  magnetic  diffusion  is  generally  much 

faster  than  thermal  diffusion  for  typical  conductors.  In  addition,  we 

have  assumed,  conservatively,  that  the  current  per  unit  height  of  rail 

in  the  vicinity  of  the  arc  is  given  by  I  /h  , 

o  a 

Let's  consider  an  arc  in  steady  state  such  that  the  ohmic 

dissipation  in  the  arc  is  just  balanced  by  the  radiation  from  its 

surface.  If  we  assume  that  the  arc  has  a  uniform  current  density,  equal 

to  I  /(wh  ),  and  a  constant  conductivity  o  and  that  it  radiates 
0  ft 

uniformly  from  its  surface,  we  can  approximate  q  by 

s 

v  w 

q  «  — - ; -  .  (5-3) 

2o  ha£a[£aw  -  £aha  +  wha] 

The  maximum  temperature  rise  at  the  rail  surface  occurs  just  as  the 

trailing  edge  of  the  arc  passes  the  point  in  question,  which  corresponds 
to  the  time 


t 


a/s 

[2a£  +  v  *] 

a  o 


v 

o 


a 


(5-4) 


where  v&  is  the  velocity  of  the  leading  edge  of  the  arc  as  it  traverses 
the  point,  and  a  is  the  acceleration  of  the  arc/projectile  system 


L  V 

2(m  +  m  ) 


(5-5) 


Combining  Eqs.  (5-1)  through  (5-5),  we  can  derive  an  expression 
for  the  velocity  below  which  we  would  expect  rail  damage  due  to  surface 
melting.  If  we  take  AT  in  Eq,  (5-1)  to  be  AT^,  the  temperature  rise 

required  to  initiate  surface  melt,  substitute  Eq.  (5-4)  for  t  and  solve 
for  vo,  we  get 


v 


o 


1_ 

2F 


( 


LI  1  II 
o  a 


m  +  m 
a  p 


(5-5) 


where 
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F  = 


np  C  k  L 
Fr  r  r 

2(m  +m  ) 
a  p 


or  h  £  (w£  +li  £  +wh  ) 
a  a  a  a  a  a 

wl 


(AT  -AT  ) 
m  R 


(5-7) 


For  copper  rails  and  parameters  typical  of  the  RM  experiment  with  an  arc 
length  of  10  cm,  Eqs,  (5-6)  and  (5-7)  yield  a  vq  on  the  order  of  2  tm/s. 

Considerable  cantion  should  be  exercised,  however,  in  applying 
Eqs.  (5-6)  and  (5-7)  to  a  given  experiment  since  the  result  is  highly 
dependent  on  the  value  assumed  for  the  arc  length.  For  instance,  for 
the  RM  experiment  F*  is  negligible  compared  to  the  first  term  in  Eq. 
(5-6).  (This  indicates  that  the  increase  in  velocity  during  the  time 
required  to  traverse  £&  is  much  less  than  vQ.)  Then,  since  m^  >>  and 

a  is  relatively  insensitive  to  £  ,  F  varies  as  £  4  for  £  >>  w  and  h  . 

a  a  a  a 

Consequently,  vq  as  predicted  by  Eq.  (5-6)  is  roughly  proportional  to 

l/£af.  Thus,  a  20- percent  increase  in  the  assumed  value  of  the  arc 

length  results  in  a  reduction  in  vq  of  nearly  a  factor  of  two.  The 

large  sensitivity  of  the  rail  temperature  rise  to  the  arc  length  also 
has  been  noted  by  Powell  [151.  In  addition,  we  note  that  vQ  is  strongly 

dependent  on  the  current.  For  a  fixed  arc  conductivity  and  arc  length, 

it  varies  as  I  4. 

o 

If  we  assume  that  the  current  is  constant  and  the  arc 
instantaneously  achieves  its  steady— state  configuration,  then  the 
temperature  rise  is  clearly  a  maximum  near  the  breech  since  the  transit 
time  there  is  a  maximum.  Thus,  one  would  expect  the  greatest  thermal 
damage  to  occur  near  the  breech.  The  velocity  vq  can  then  be 

interpreted  as  the  velocity  with  which  the  projectile  must  be  injected 
in  order  to  prevent  melting  of  the  rails  at  the  breech. 

Experimental  data  does  indeed  suggest  that  the  damage  to  the 
rails  in  arc-driven  rail  guns  is  greatest  near  the  breech.  It  should  be 
mentioned,  however,  that  the  steady-state  assumption  is  probably  not 
valid  near  the  breech  [5).  Generation  of  an  arc  by  the  exploding-wire 
technique  leads  initially  to  a  highly  non-steady  situation,  with  the 
plasma  internal  pressure  causing  the  arc  to  expand  against  the  Lorentz 
force  which  acts  to  force  the  arc  against  the  projectile.  The  details 
of  this  interaction,  which  has  yet  to  be  examined  either  experimentally 
or  analytically,  will  surely  influence  the  energy  transport  to  the  rails 
near  the  breech. 
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VI.  SUMMARY  AND  DISCUSSION 


In  this  report,  we  have  developed  an  equation  for  the  effective 
inductance  per  unit  length  for  rail  guns  with  unequal  arc  and  rail 
heights.  The  equation  provides  a  simple  technique  for  estimating  the 
net  accelerating  force,  based  solely  on  the  current  and  the  gun 
geometry. 

In  addition,  we  have  developed  a  technique  for  modifying  previous 
models  of  arc  dynamics  in  rail  guns  to  account  for  the  finite  height  of 
the  rails.  For  typical  rail  gun  geometries,  the  model  suggests  that  the 
pressure  profile  for  f inite-height  rails  can  be  well  approximated  by 
merely  scaling  the  pressure  profile  for  inf inite-height  rails,  the 
scaling  factor  being  a  function  of  the  effective  inductance  for  the  two 
geometries.  Exploiting  this  result  in  a  simple  arc  model  for  the  RM 
experiment,  we  find  that  use  of  the  finite-height  model  leads  to 
considerably  lower  predictions  for  the  arc  pressure,  density, 
acceleration,  and  arc  mass  when  compared  to  the  predictions  for  the 
infinite  height  model  at  the  same  arc  length.  On  the  other  hand,  the 
use  of  the  finite-height  rail  gun  model  leads  to  a  higher  prediction  for 
the  arc  muzzle  voltage.  The  trends  indicated  by  the  infinite-height 
model,  however,  remain  unchanged  when  the  effect  of  finite-height  rails 
is  included. 

In  future  work,  the  finite-height  model  developed  in  this  report 
should  be  incorporated  into  the  more  detailed  one-dimensional  model 
described  in  Ref.  5.  This  should  provide  more  accurate  predictions  of 
the  properties  of  arcs  in  arc-driven  rail  guns.  Furthermore,  the 
properties  predicted  by  the  finite-height  model  should  he  used  in  any 
future  studies  of  arc  stability. 

A  review  of  the  recent  literature  has  revealed  that  models  have 
been  proposed  for  analyzing  mechanical  and  thermal  damage  to  the  rails 
in  arc-driven  rail  guns.  Although  the  assumptions  used  iu  the  models 
appear  reasonable,  there  is  insufficient  experimental  evidence  to 
determine  their  validity.  The  thermal  damage  model,  in  particular,  is 
sensitive  to  parameters,  such  as  the  arc  length,  which  are  difficult  to 
estimate  accurately.  Although  this  limits  its  utility  for  reproducing 
experimental  data,  the  model  is  valuable  nonetheless  for  indicating  the 
measures  that  can  be  taken  to  reduce  thermal  energy  transfer  to  the 
rails.  For  instance,  the  model  predicts  that  increasing  the  arc  length 
(by  increasing  the  arc  mass,  for  example)  will  lead  to  considerably 
lower  beat  transfer  to  the  rails. 

One  difficulty  with  the  thermal  model  is  the  assumption  that  the 
arc  is  characterized  by  a  steady  state.  This  assumption  is  probably  not 
justified  near  the  breech  where  thermal  damage  to  the  rails  is  greatest. 
Consequently,  future  studies  of  rail  damage  due  to  the  arc  should 
account  for  the  time  dependent  processes  occurring  during  the  generation 
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of  the  arc.  In  addition,  the  analysis  has  been  conducted  only  for  a 
single-shot  rail  gun.  A  relatively  straightforward  extension  of  the 
model  to  analyze  the  cumulative  effects  of  several  shots  would  be 
particularly  useful  for  determining  how  rail  heating  will  influence  the 
operation  of  a  repetitively  fired  gun. 

Finally,  experimental  work  needs  to  be  done  to  characterize  the 
arc  and  to  provide  estimates  for  the  extent  and  type  of  damage  to  the 
rails.  Some  initial  work  in  this  area  has  been  reported  [16,  17]  and 
these  efforts  hopefully  will  be  continued.  The  work  of  Jamison  et  al  at 
BRL  is  particularly  noteworthy.  On  the  basis  of  tracks  left  on  the 
conducting  rails  and  light  emission  from  the  arc,  they  conclude  that 
fine  structure  exists  in  the  current  distribution  in  the  arc.  The 
existence  of  preferred  current  paths  and  local  *hot  spots'  in  the  arc 
would  significantly  affect  heat  transfer  to  the  rails.  Clearly, 
additional  study  of  this  phenomenon,  both  experimentally  and 
theoretically,  would  be  desirable. 
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APPENDIX  A 


In  this  Appendix*  we  show  that  the  net  electromagnetic  force  on 
the  arc/projectile  system  is  independent  of  the  current  distribution  in 
the  arc  for  the  f inite-height-rail  model  of  Sec.  III* 

The  electromagnetic  force  is 


F  =  wh  8  f  <J  >  <B  >  d£  . 
m  a  J0  y  z 

Substituting  (3-21)  into  Eq.  (A-l)  and  using  the  definitions 


yields 


where 
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(A-l) 


(A-2) 


(A-3) 


(A-4) 


If  we  integrate  R  by  parts  with  respect  to  £*  then  repeat  the  procedure 
with  respect  to  and  add  the  results*  we  obtain 

a 

R  =  |d(1)  J  du/d$  |  Q[ (£-1)  Ca/w]  +  Q((l-§)  £ft/w]  |  d$  .  (A-5) 

Substituting  Eq.  (3-23)  into  Eq.  (A-5)  and  integrating  yields 


R  =  [S(«>  -  S(0)J  u(l) 


(A-6) 


Replacing  R  by  this  result  in  Eq.  (A-3)  and  utilizing  the  definition  of 
u(l)  given  in  Eq.  (A-2)  yields  our  final  result* 
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<A-7) 


pwh  j  »  r  s<«o  -  s(o) 

Fm  *  2  L1  ~  nV 

Thus,  the  net  electromagnetic  force  depends  only  on  the  current  per  unit 
height  of  rail  and  the  rail  gun  geometry. 


50 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 


12  Administrator  1 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22314 

1  Office  Under  Secretary  of  Defense 

Research  5  Engineering  6 

ATTN:  Mr.  Ray  Thorkildsen 
Room  3D 10 89,  The  Pentagon 
Washington,  DC  20301 

1  Deputy  Under  Secretary  of  Defense 
Research  §  Engineering 
Room  3E114,  The  Pentagon 
Washington,  DC  20301 

1? 

3  Director 

Defense  Advanced  Research  Projects 
Agency 

ATTN:  Dr.  Joseph  Mangano 

Dr.  Gordon  P.  Sigman 
Dr.  Harry  Fair  1 

1400  Wilson  Boulevard 
Arlington,  VA  22209 

1  Office  of  Assistant  Secretary  of 
the  Army 

ATTN:  RDA,  Dr.  Joseph  Yang 

Room  2E672,  The  Pentagon 

Washington,  DC  20310  1 

1  HQDA  (DAMA-ARZ-A/Dr.  Richard  Lewis] 
Washington,  DC  20310 

1  Commander 

US  Army  Materiel  Development  § 

Readiness  Command  1 

ATTN:  DRCDMD-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

1  Commander 

US  Army  Materiel  Development  §  1 

Readiness  Command 
ATTN:  DRCLDC,  Mr.  Langworthy 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333 


Commander 

US  Army  Armament,  Munitions  & 

Chemical  Command 

Armament  Research  §  Development  Center 
ATTN:  DRSMC-TDC (D) 

Dover,  NJ  07801 

Commander 

US  Army  Armament,  Munitions  § 

Chemical  Command 
ATTN:  DRSMC-TSS (D) 

DRSMC-LCA(D) ,  Mr.  J.A.  Bennett 
Dr.  T.  Gora,  Dr.  P.  Kemmy 
DRSMC-SCA(D),  Mr.  W.R.  Goldstein 
Dover,  NJ  07801 

Commander 

US  Army  Armament,  Munitions  § 

Chemical  Command 
ATTN:  DRSMC-LEP-L (R) 

Rock  Island,  IL  61299 

Director 

Benet  Weapons  Laboratory 

Armament  R&D  Center 

US  Army  Armament ,  Munitions  § 

Chemical  Command 
ATTN :  DRSMC -LCB -TL 
Watervliet,  NY  12189 

Commander 

US  Army  Aviation  Research  § 

Development  Command 
ATTN :  DRDAV-E 
4300  Goodfellow  Boulevard 
St.  Louis,  MO  63120 

Director 

US  Army  Air  Mobility  Research  $ 
Development  Command 
Ames  Research  Center 
Moffett  Field,  CA  94035 

Commander 

US  Army  Communications  Research  6 
Development  Command 
ATTN:  DRSEL-ATDD 
Fort  Monmouth,  NJ  07703 


51 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 
1  Commander 

US  Army  Electronics  R§D  Command 
Technical  Support  Activity 
ATTN:  DELSD-L 
Fort  Monmouth,  NJ  07703 

1  Commander 

US  Army  Missile  Command 

ATTN:  DRSMI-R 

Redstone  Arsenal,  AL  35898 

1  Commander 

US  Army  Missile  Command 
ATTN:  DRSMI-YDL 
Redstone  Arsenal ,  AL  35898 

1  Commander 

US  Army  Tank  Automotive  Command 
ATTN :  DRSTA-TSL 
Warren,  MI  48090 

1  Director 

US  Army  TRADOC  Systems  Analysis 
Activity 
ATTN:  ATAA-SL 
White  Sands  Missile  Range 
NM  88002 

2  Commandant 

US  Army  Infantry  School 
Am:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905 

2  Commander 

US  Army  Research  Office 
ATTN:  Dr.  Fred  Schmiedeshoff 
Dr.  M.  Ciftan 
PO  Box  12211 

Research  Triangle  Park,  NC  27709 
1  Commander 

Naval  Air  Systems  Command 
ATTN:  John  A.  Reif,  AIR  350B 
Washington ,  DC  20360 


No.  of 

Copies  Organization 
4  Commander 

Naval  Surface  Weapons  Center 
ATTN:  Mr,  H.B.  Odom,  Code  F-12 
Dr,  M.F.  Rose,  Code  F-04 
Mr.  P.T.  Adams  and 
Mr,  D.L.  Brunson,  Code  G-35 
Dahlgren,  VA  22448 

2  Commander 

US  Naval  Research  Laboratory 
ATTN:  Dick  Ford,  Code  4774 
Dr.  I.  Vitkovitsky 
Washington,  DC  20375 

1  HQ  AFSC/SD0A 

ATTN:  XRB  (CPT  Dennis  Kirlin) 

Andrews  AFB,  MD  20334 

2  AFATL/DLDB  (Lanny  Burdge,  Bill  Lucas) 
Eglin  AFB,  FL  32542 

1  AFATL  (Richard  Walley) 

Eglin  AFB,  FL  32542 

1  AFWL  (Dr.  William  L.  Baker) 

Kirtland  AFB,  NM  87117 

1  AFWL/NTYP  (John  Generosa) 

Kirtland  AFB,  NM  87117 

1  AFWL/ SUL 

Kirtland  AFB,  NM  87117 

1  AFAPL  (Dr.  ChaTles  E.  Oberly) 

Wright- Patterson  AFB,  OH  45433 

1  AF AWL/POOS- 2  (CPT  Jerry  Clark) 

Wright -Patterson  AFB,  OH  45433 

1  Director 

Brookhaven  National  Laboratory 
ATTN:  Dr.  James  R.  Powell,  Bldg  129 
25  Brookhaven  Avenue 
Upton,  NY  11973 

2  Director 

Lawrence  Livermore  National  Laboratory 
ATTN:  Dr.  R.S.  Hawke 
L-156,  Box  808 
Livermore,  CA  94550 


52 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 


No.  of 

Copies  Organization 

2  Director 

Los  Alamos  Scientific  Lab 
ATTN:  Dr.  Clarence  M.  Fowler,  MS970 
Dr.  Denis  R.  Peterson,  MS985 
Los  Alamos,  NM  87545 

3  NASA-LeRC-MS  501-7 
ATTN:  Bill  Kerslake 

Frank  Terdan 
Mike  Brasher 
21000  Brookpark  Road 
Cleveland,  OH  44135 

1  Boeing  Aerospace  Company 
ATTN:  J.E.  Shrader 
PO  Box  3707 
Seattle,  WA  94124 

1 •  General  Dynamics 

ATTN:  Dr.  Jaime  Cuadros 
P0  Box  2507 
Pamona,  CA  91766 

3  GT  Devices 

ATTN:  Dr,  Derik  Tidraan 

Dr.  Shyke  Goldstein 
Dr.  Neils  Winsor 
570 5-A  General  Washington  Drive 
Alexandria,  VA  22312 

1  Physics  International 
ATTN:  Dr.  A.L.  Brooks 
2700  Merced  Street 
San  Leandro,  CA  94577 

2  R$D  Associates 
ATTN:  Mr.  Ronald  Cunningham 

Dr.  Peter  Turchi 
PO  Box  9695 

Marina  del  Rey,  CA  90291 

1  Science  Applications,  Inc. 

ATTN:  Dr.  Jad  H.  Batteh 
1503  Johnson  Ferry  RD,  Suite  100 
Marietta,  GA  30062 


Science  Applications,  Inc. 

Corporate  Headquarters 
ATTN:  Dr.  Frank  Chilton 
1250  Prospect  Plaza 
La  Jolla,  CA  92037 

1  Unidynamics/Phoenix,  Inc. 

ATTN:  W.R.  Richardson 
P0  Box  2990 
Phoenix,  AZ  85062 

2  Vought  Corporation 
ATTN:  William  B.  Freeman 

Charles  Haight 
P0  Box  225907 
Dallas,  TX  75265 

2  Westinghouse  Research  8 
Development  Laboratory 
ATTN:  Dr.  Ian  R.  McNab 
Dr.  Y.  Thio 
Gateway  3 

Pittsburgh,  PA  15222 

2  Massachusetts  Institute  of  Technology 
Francis  Bitter  National  Magnet  Lab 
ATTN:  Dr.  Henry  H.  Kolm, 

Dr.  Peter  Mongeau 
170  Albany  Street 
Cambridge,  MA  02139 

1 •  University  of  Texas 

Center  of  Electroraechanics 
ATTN:  Mr.  William  F.  Weldon 
167  Taylor  Hall 
Austin,  TX  78712 

Aberdeen  Proving  Ground 

Dir,  USAMSAA 
ATTN:  DRXSY-D 

DRXSY-MP,  H.  Cohen 
Cdr ,  USATECOM 
ATTN:  DRSTE-T0-F 

Cdr,  CRDC,  AMCCOM,  Bide  E3516 
ATTN:  DRSMC-CLB-PA 
DRSMC-CLN 
DRSMC-CLJ-L 


53 


THIS  PAGE  WAS  INTENTIONALLY  LEFT  BLANK 


USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/ contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


- -  -  -  -  - - - FOLD  HERE  - 

Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRSMC-BLA-S  (A) 

Aberdeen  Proving  Ground,  MD  21005 


OFFICIAL.  BUSINESS 

PENALTY  FOR  PRIVATE  USE.  S3O0 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  12062  WASHINGTON, 0C 
POSTAGE  WILL  BE  PAIO  BY  DEPARTMENT  OF  THE  ARMY 


Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DR SMC-BLA-S  (A) 

Aberdeen  Proving  Ground,  MD  21005 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


